clc; clear; close all

% 参数设置
b2 = 2*pi;
c  = 11.5;
a_list = linspace(0.33, 5, 400);  % 参数a区间
N = 800;          % 总步数（含暂态）
transient = 300;  % 暂态步数
delta = 1e-7;     % 差分步长  
    LEs = LE(a, b2, c, x0, y0, z0, N, transient, delta);
    LEs_all(ai,:) = sort(real(LEs), 'descend');
end

% 绘图
figure
plot(a_list, LEs_all(:,1), 'r', 'LineWidth', 1.2); hold on
plot(a_list, LEs_all(:,2), 'g', 'LineWidth', 1.2);
plot(a_list, LEs_all(:,3), 'b', 'LineWidth', 1.2);
xlabel('a'); ylabel('LEs');
legend('LE_1', 'LE_2', 'LE_3')
